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Abstract 

We revisit two old and apparently little known papers by Basuev m si and show that the 
results contained there yield strong improvements on current lower bounds of the convergence 
radius of the Mayer series for continuous particle systems interacting via a very large class 
of stable and tempered potentials which includes the Lennard-Jones type potentials. In 
particular we analyze the case of the classical Lennard-Jones gas under the light of the 
Basuev scheme and, using also some new results [3B] on this model recently obtained by one 
of us, we provide a new lower bound for the Mayer series convergence radius of the classical 
Lennard-Jones gas which improves by a factor of the order 10 5 on the current best lower 
bound recently obtained in [T9] . 


1 Introduction 

The possibility to obtain the equation of state for a non-ideal gas only from first principles, i.e. once 
given the microscopic interaction between its molecules, has been a subject of intense investigation 
by the mathematical physics community and several rigorous results on this issue has been obtained, 
especially in the decade of the sixties, concerning mainly systems of classical particles interacting 
via a pair potential. 

In this context, Mayer, [zn [22 ] and Mayer and Mayer [23] . in their seminal works in the forties, 
were able to obtain an explicit representation of the pressure of a non ideal gas whose particles 
interact via a pair potential. Such an expression, known nowadays by the name of Mayer series, is 
a formal power series in terms of the activity (a quantity trivially related to the chemical potential) 
whose coefficients (called the Ursell coefficients) can be computed once given the microscopic pair 
potential. 

The Mayer series represents an explicit expression of the equation of state of a non ideal gas 
deduced from first principles (i.e. only from the microscopic pair interaction between particles). 
However the series was at that time only formal, in the sense that nothing could be said about 
its convergence due to the intricate combinatorial structure of the Ursell coefficients. Even so, 
the Mayer series (truncated at some order) has been successfully used since then by chemists 
and physicist to extract useful information about the behavior of real gases. This was a strong 
motivation to prove its convergence, at least in the region of low density or high temperature, where 
every classical system of particles is supposed to be in the gas state. Such an achievement would 
provide a firm ground to the interpretation of the Mayer series as the rigorous and exact equation 
of state for any non-ideal gas. 
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It took more than twenty years after Mayer’s results to obtain the first rigorous results on the 
convergence of the Mayer series. The first breakthrough, due to Groeneveld m, was a proof of the 
convergence of the Mayer series in the low-density/high-temperature phase for gases of particles 
that interact via a purely repulsive pair potential. The latter was a quite dramatic restriction, 
which ruled out most, if not all, pair potentials modeling realistic continuous particle systems in 
physics. However, just one year later, Penrose [25, 26] and independently Ruelle [34J [35] obtained 
the astonishing result that the Mayer series of continuous systems of particles was actually an 
analytic function of the temperature and activity in the low density/high temperature phase. The 
result holds for an enormous class of pair interactions (the so-called stable and tempered pair 
potentials, see ahead for the definitions), which included practically all known examples of pair 
potentials for realistic gases. More than this, Penrose and Ruelle were able to provide a lower 
bound for the convergence radius of the Mayer series, which stands among the best available in the 
literature till nowadays. 

The mathematical methods developed by Penrose and Ruelle to get their results were based on the 
so called Kirkwood-Salzburg equations (see e.g. Sec. 4.2.1 in [33] and references therein). Such 
techniques do not face directly the structure of the general n th order Ursell coefficient of the Mayer 
series (a well behaved upper bound for it, e.g., such as C n for some constant C, would provide 
immediately the analyticity of the series for low densities). Rather, these methods are based on the 
analysis of an infinite set of relations relating the n-point correlation function of the system (which 
can also be expressed as a series similar to the Mayer series) to the m-point correlation functions 
with m < n). 

Despite the undoubting success of the methods based on Kirkwood-Salzburg equations, the appar¬ 
ently great difficulty to obtain a direct bound for the Mayer coefficients starting from their explicit 
expression in terms of sum over connected graphs remained an intriguing open question. 

The first result in this direction was obtained by Penrose in 1967 m who was able to obtain 
bounds on the Mayer coefficients, as powerful as those of Penrose-Ruelle in 1963 ['26] 134) . directly 
via a resummation of the connected graphs appearing in the Ursell coefficients in terms of trees: it 
was the first example of tree graph identity. This very nice result, however, was restricted to pair 
potentials which had not only to be stable and tempered but also possess a hard-core. The Penrose 
tree-graph identity, somehow forgotten for thirty years, was recently rescued in mm and [HJ 
where it has been utilized to improve the convergence region of the abstract polymer gas and the 
gas of hard spheres in the continuum respectively. 

In 1978, Brydges and Federbush [5] developed a new tree-graph identity (different from that of 
Penrose) through which they were able to obtain bounds stronger than those of Penrose-Ruelle 
for a subclass of the stable and tempered potentials. Pair potentials in this subclass needed to be 
absolutely summable. This was a quite severe restriction since it left out the important class of the 
Lennard-Jones type potentials. However a quite popular potential among chemists was included: 
the so called Morse Potential, and in this case the improvement was really consistent (see e.g. [24]). 
The original identity presented in [5] was successively reworked and extended in several papers (H. 
0 , 0 , m, 0 ) and its quite robust structure has been successfully used mainly for applications in 
constructive field theory, especially during the eighties (see e.g. m and references therein). 

There have been also more recent developments of the Brydges-Federbush tree-graph identity which 
are directly related with the present paper (see e.g. eu, m and [29]). In particular, via the 
Brydges-Federbush identity, Morais et al. [24] were able to improve the Penrose-Ruelle region of 
analyticity of the Mayer series for (some particular cases of ) Lennard-Jones type potentials. Based 
on the ideas of [24] (and using also some results in [20]) de Lima and Procacci [19] . improved 
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strongly the convergence radius of the classical Lennard-Jones potential. 

In 1978, the same year of the original paper by Brydges and Federbush, A. G. Basuev published (in 
the Soviet Journal Theoreticheskaya i Mathematicheskaya Fizika) a paper [3] providing a stability 
criterion for a quite large class of pair potentials (which included the Lennard-Jones type potentials). 
If such criterion holds then the potential can be split into a short distance non-negative compactly 
supported term (the repulsive part) plus an absolutely integrable term (incorporating the possible 
attractive long range tail of the potential). The stunning property of the absolutely integrable 
part is that it has the same stability constant as the full potential. This very interesting result 
was used by the same author one year later [3] to obtain bounds for the Mayer coefficients for 
this class of potentials (which, we repeat, is very large and possibly includes nearly all reasonably 
physical example). These bounds appeared to be strongly better than those obtained by Penrose 
and Ruelle. Even more surprising, to get these bounds the author constructed an original tree graph 
identity, different from that of Penrose of 1967 m as well as that of Brydges-Federbush of 1978 
[5]. Apparently the author was not aware at the time of the paper [5], so he did not try to compare 
his result with that of Brydges and Federbush. If he did, he would arrive at the conclusion that his 
methods yield (in the worst of the hypothesis) the same improvements given in [5] for absolutely 
summable pair potentials but of course are able to provide fantastic bounds for non absolutely 
summable pair potential which were at the time out of reach with the Brydges-Federbush methods. 
These two works by Basuev, which should be cited in any paper dealing with the convergence of 
the Mayer series for continuous systems, were instead nearly completely overlooked. Actually the 
Basuev stability criterion presented in [3] was cited in some detail once [;2Ei|. We were not able 
to find any reference at all in the literature to the Basuev tree-graph identity and the consequent 
improvements on the Mayer coefficient bounds presented in the second paper |3J. 

In the present work we revisit these two papers by Basuev with the intent to rescue the important 
results there contained, which apparently passed completely unobserved. In particular, we show 
how these results yield improvements on all known results for a wide class of pair potentials, 
focusing in particular on the important and physically relevant subclass of the Lennard Jones type 
potentials. We also provide new results for the specific case of the classical Lennard-Jones potential 
V(r) = 1/r 12 — 2/r 6 which strongly improve on the results given by two of us in a recent paper |19]. 
To get this last result on the specific case of the Lennard-Jones potential, we used the very recent 
results obtained by one of us [38] on the stability constant and minimal inter-particle distance in 
lowest energy configurations of the Lennard-Jones gas. 

The rest of the paper is organized as follows. In section 2 we introduce notations and model. We 
recall the definition of Mayer series and Mayer coefficients. We introduce the concepts of stability 
and temperedness, and the class of potentials satisfying the Basuev criterion (Definition [3]). We 
finally recall the main results on Mayer series convergence by Penrose and Ruelle (Theorem [T]), and, 
the new recent results by Morais, Procacci and Scoppola (Theorem[2]). In section 3 we (re)introduce 
the Basuev tree graph identity (Theorem [3]) which we (re)prove in Appendix A. In section 4 we 
(re)state the Basuev criterion (Theorem[3|), whose proof is given in Appendix B, and (re)derive from 
the Basuev tree graph identity the bounds on the Ursell coefficient for particle systems interacting 
via a Basuev pair potential. Finally, in Section 5 we present the estimates of the Mayer coefficients 
and convergence radius for Lennard-Jones type potentials and classical Lennard-Jones potential, 
showing how this estimates improve on the recent bounds obtained in |24] and [19] . 
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2 The model 


We consider a continuous system of particles in the d-dimensional Euclidean space W l . Denote by 
X{ € the position vector of the z th particle of the system and by \xi\ its Euclidean norm. We 
suppose hereafter that particles interact through a distance dependent pair potential V(r) with 
r € [0,+ 00 ). Given a configuration (x\,... ,x n ) £ M dn of the system such that n particles are 
present, the energy U{x 1 ,... ,x n ) of this configuration is defined as 

U(xi,...,x n ) = ^ V{\xi-Xj |) 


The statistical mechanics of the system is governed by the following partition function in the Grand 
Canonical Ensemble 


00 \n r r 

S a (A,/3) = V'— / dx\ ... / d,x n e~P u{xi '-' Xn) 

n! J A J A 


where A C is the region where the system is confined (which hereafter we assume to be a 
d-dimensional cube with sides of length L and denote by |A| its volume), A is the activity of the 
system and /3 is the inverse temperature. 

Thermodynamics is recovered by taking the logarithm of the partition function. Namely, the 
pressure P\ and the density p\ of the system at fixed values of the thermodynamic parameters 
inverse temperature /3, fugacity A and volume A, are given respectively by the following formulas 

PA = M logHA(A,/3) (2 - 2) 

PA = tkt Jx log Ha ^ A ’ ^ ( 2 - 3 ) 

As we will see below, both P\ and p\ are given in terms of a power series of the fugacity A (with 
coefficients depending on j3). Inverting (2.3) by writing A as a power series of the density p and 
plugging this series into (2.2) one obtains the so called Viral series, i.e. the pressure as a function 
of the density and temperature, or, in other words, the equation of state of a gas whose particles 
interact via the pair potential V. 

The dependence of the pressure Pa and density p\ on the volume A should be a residual one. In 
fact, one may think to increase the volume A of the system keeping fixed the value of the fugacity 
A and the inverse temperature (3. It is then expected that pressure and density of the system do 
not vary significantly. It is usual to think that the box A can be made arbitrarily large (which is 
the rigorous formalization of macroscopically large) and define 


P = A 1 TooM 10gEA(A ’ /3> (2A) 

p = A 1 i5ojidJd og = A(A ' ,3) (2 - 5) 

The limit A —>• 00 (here A —>• 00 means that the side L of the cubic box A goes to infinity) is usually 
called by physicists the thermodynamic limit and the exact thermodynamic behavior of the system 
is in principle recovered at the thermodynamic limit. 
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It is thus crucial for the study of the thermodynamic behavior of the system to have available an 
explicit expression for |A| _1 logEA(A, f3) and to be able to control its behavior as A —>• oo. 

A very well known and old result due to Mayer and Mayer [23j (see there chapter 13, p. 277-284) 
states that the quantity logHA(/3, A), with Ea(/3, A) being the partition function defined as in (2.1), 
can be written in terms of a formal series in powers of A. Namely, 


In “a(A,/3) 


|A| 


A + ^C n (/3,A)A ri 

n =2 


( 2 . 6 ) 


where 


C n (P, A) 



e ~PV{\xi-Xj\)\ _ i 


(2.7) 


with G n denoting the set of all connected graphs with vertex set [n] = {1,2 ,n} and E g denoting 
the edge set of g E G n . 

The power series (2.6) is the famous Mayer series, the factors C n (j3, A) are usually called Mayer 
(or Ursell) coefficients and the factors <f)p(xi,... ,x n ) are usually called the Ursell functions (see 
e.g. [33] and references therein). 

In order for the limits (2.4) and (2.5) to exist, the dependence on the volume |A| of the Mayer 
coefficients C n (f3, A) given by (2.7) must be only a residual one. Indeed, via the expression (2.7) it 
is not difficult to show that C n (/3, A) is bounded uniformly in A and, for every n E N, the limit 


C n (P)= lim C n (/3, A) 

|A|-kx> 


( 2 . 8 ) 


exists and it is a finite constant (see e.g [33]). On the other hand it is a very difficult task to 
obtain an upper bound on \C n (/3, A)| which is at the same time uniform in the volume |A| and 
well-behaved in n : i.e. an upper bound of the form 

\C n {p,A)\<(Cp) n (2.9) 

with C /3 some positive constant. Such an estimate would immediately yield a lower bound for the 
convergence radius of the Mayer series of the pressure uniform in the volume. The difficulty to 
obtain a bound of type (2.9) directly from the expression (2.7) of the Mayer coefficients stems from 
the fact that the cardinality of the set G n (the connected graphs with vertex set [n]) behaves very 
badly with n (it is no less than 2("= 1 )( n_2 )/ 2 , which is the number of connected graphs containing 
a fixed tree). Therefore, it seems more than reasonable to impose some conditions on the potential 
V and at the same time it seems inevitable to look for some hidden cancelations in the expression 
(2.7). 

During the sixties it became clear that the minimal conditions one has to impose on the potential V 
to get bounds for the Ursell coefficients of the type (2.9) were basically stability and temperedness, 
whose definitions we recall below. 


Definition 1 (Stability) A pair potential V is said to be stable if there exists a finite C > 0 such 
that, for all n & N and all (x \,..., x n ) E M. dn , 

U(xi,...,x n )= ^ V(\xi - Xj\) >-nC (2.10) 

l<i<j<n 
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Definition 2 (Temperedness) A pair potential V is said to be tempered if there exists rg > 0 
such that 


'\x\>r 0 


dx\V{\x\)\ < oo 


( 2 . 11 ) 


Stability and temperedness are actually deeply interconnected and the lack of one of them always 
produces non thermodynamic or catastrophic behaviors (see e.g. the nice and clear exposition 
about the role of stability and temperedness in reference m)- As far as stability is concerned, 
observe that the grand-canonical partition function defined in (2.1) is a holomorphic function of A if 
the potential V is stable. Moreover, under very mild additional conditions on the potential (upper- 
semicontinuity) it can be proved that the converse is also true (see Proposition 3.2.2 in [33]). In 
other words /3) converges if and only if the potential V is stable. So, in some sense, stability 

is a condition sine qua non to construct a consistent statistical mechanics for continuous neutral 
particle system^]. Concerning temperedness, while it is reasonable to expect that the potential 
decays to zero at large distances, on the other side it is not completely evident to see why the 
potential should decay in an absolutely integrable way. E.g., there are examples of charged particle 
systems, such as the dipole gas, which interact through a non-tempered pair potential and have a 
convergent Mayer series. We denote 

B n = sup -~U(xx n ) (2.12) 

(xi,...,Xn)(zM. dn 71 

B = supB n (2.13) 

n> 2 

We also define, for later use 

_ fl _ _ 

B n = - -B n B = supB n (2.14) 

n - 1 n > 2 


As will see below, the quantities B and B play a crucial role in the estimates of the convergence 
radius of the Mayer series. Note that B , when finite, is actually the best constant C in (2.10). 
Note also that temperedness of V implies that B is non-negative and B = 0 if and only if V > 0 
(i.e. “repulsive” potential). The non-negative number B defined in (2.13), when finite, is known as 
the stability constant of the potential V. The slight variant B defined in (2.14) has been originally 
introduced in [3j. Of course, by definition B > B and B and B, if not coinciding, should be 
very close for most of the reasonable stable pair potentials. As noted by Basuev [3J, it holds that 
B < B < j^B for any stable and tempered three-dimensional pair potential which is definitively 
negative at sufficiently large distances. 

In the present paper we work with a proper subclass of the tempered and stable pair potentials 
which was first proposed by Basuev in [3]. This class is sufficiently large to embrace practically 
every known example of physically relevant pair potential and we will see in Section [4] that a 
potential in this class has a key property (Theorem [4] ahead) which is the main ingredient toward a 
very effective bound of the form (2.9) for the Ursell coefficient (2.7). The class is defined as follows. 


1 For systems of charged particles stability is no longer a necessary condition, see e.g. the example of the 
two-dimensional Yukawa gas in [7] or [8]. 
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Definition 3 A tempered pair potential V is called Basuev if there exists a > 0 such that 


V{r) > V(a) > 0 


for all r < a 


and 

with 


V(a) > 2/r(o) 


/i(o) = sup (|xi|) 

n€N, (xi ,...,x n )eR dri 
\ x i~Xj\>a V{i,j}eE n 


(2.15) 

(2.16) 

(2.17) 


where V is the negative part of the potential, i.e. V = \ |V"| — V 

As mentioned in the introduction, the best rigorous bound on |C n (/3, A)| available in the literature 
for stable and tempered pair potentials is that given by Penrose and Ruelle in 1963 


Theorem 1 (Penrose-Ruelle) Let V be a stable and tempered pair potential with stability con¬ 
stant B. Then the n-order Mayer coefficient C n (/3, A) defined in (2.7) is bounded by 

\C n {P,A)\ < e WB(n-l) n n-2 [C(P)] n 1 ( 2 . 18 ) 

n\ 

where 

C(/3) = [ dx le-^l*!) - 1| 

J R d 

Therefore the Mayer series (2.6) converges absolutely, uniformly in A for any complex X inside the 
disk 

l A l < e 2f)B+l C (p) ( 2 ‘ 19 ) 

The Penrose-Ruelle bound has been improved for some restricted classes of stable and tempered 
pair potentials. Brydges and Federbush [5] gave an improvement of the Penrose-Ruelle bound as 
far as absolutely summable pair potentials are considered and it can be seen that such improvement 
is quite strong for several physically relevant pair potentials (see e.g. the example of the Morse 
potential treated in [24]). Recently, via an extension of the Brydges-Federbush identity, Morais et 
al. [21] improved the Penrose-Ruelle bound, by considering pair potentials that can be written as 
a sum of a positive part plus an absolutely integrable stable part. 


Theorem 2 (Morais-Procacci-Scoppola) Let V be a stable and tempered pair potential with 
stability constant B. Let V be such that V = <hi + < I ) 2 , with <3?i > 0 and <f >2 absolutely integrable and 
stable with stability constant B. Then the n-th order Mayer coefficient C n (/3, A) defined in (2.7) 
admits the bound 

\C n (/3, A)| < e^ n n n ~ 2MMII (2.20) 

n\ 


where 



1| +/3|*h 2 (M)| 


( 2 . 21 ) 
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Consequently, the Mayer series (2.6) converges absolutely for all complex activities A such that 


|A|< 


1 

eh^+ 1 C(/3) 


( 2 . 22 ) 


Remark. We draw the attention of the reader to the fact that in expression (2.22) for the lower 
bound of the convergence radius there appears the stability constant B of the absolutely summable 
part <f >2 of the potential in place of the stability constant B of the full potential V as in (2.19). 
This is somehow unpleasant since in general B > B, so that bound (2.22) has chance to be an 
improvement of the Penrose-Ruelle bound (2.19) only if B < 2 B. Actually in [24] the authors did 
produce some examples of pair potentials that can be decomposed as above, with B = B. Moreover 
in m it is shown that also the classical Lennard-Jones potential (see ahead for the definition) has 
this property. As we will, all these examples fall in the class of the Basuev potentials according to 
Definition [3] and are thus covered by Theorem [5] in section 4 below. 


3 The Basuev tree-graph identity 

In this section we present the Basuev tree-graph identity. As the Penrose identity (formula (6) in 
m) and the Brydges-Federbush identity (formula (50) in [5]), it is an alternative expression of the 
Ursell coefficients (2.7) in terms of a sum over trees rather than connected graphs. As mentioned 
in the introduction and as we will see ahead with all details, this identity permits to get rid of 
the combinatorial problem by directly bounding the Ursell coefficients starting from their explicit 
expressions (2.7). The advantage of the Basuev tree-graph identity given in [4] respect to the 
previous ones given by Penrose m and by Brydges and Federbush [5] is twofold: it permits to 
treat a much larger class of pair potentials and it permits to obtain better bounds. It is worth to 
say, that the Basuev identity can be deduced as a particular case of most recent versions mmn ° f 
the Brydges-Federbush tree-graph identity, e. g. such as that presented in Theorem VIII.3 of [8]. 
The Basuev formula is essentially algebraic and combinatorial. In order to introduce it, we need 
to give some preliminary notations. 

Given a set X, we denote by P(A) the set of subsets of X (i.e. P(A^) is the power set of X) and 
by P*(A) the set of non-empty subsets of X (i.e. P*(A) = P(A) \ {0}). 

Set partitions. Given a finite set X , we denote by 11(A) the set of all partitions of X. Namely, 
an element n € n(A) is, for k = 1,... |A|, a collection it = {m, ... , a*,} of non-empty pairwise 
disjoint subsets of X such that U= X. The elements Oj {i = 1,..., k) of the collection forming 
7r are called the blocks of 7r. Given a partition ir = {ai,... , a*.} of X we set 17r| = k. For fixed 
k = 1,..., |A| we denote by II/ C (A) the set of all partitions of X into exactly k blocks, or, in other 
words, Ilfc(A) is the set of all ir € n(A) such that 17r| = k. 

Let n € N. We remind the notation [ n ] = {1,2,... ,n}. We also denote by E n the set of all 
unordered pairs in [n]. 

Definition 4 A pair interaction in the set [n] is a map V : E n —» MU{+oo} that associates to any 
unordered pair in {i,j} € E n a number Vij = V({i, j}) = Vji with values in M U {+oo}. 
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Definition 5 Given a pair interaction V in [to], the energy induced by V is a map U : P([to]) —> M 
defined as follows. For any subset X of [to] 


u(x) = ]T Vij 

{i,j} CX 


(3.1) 


Definition 6 Let (5 € (0, +oo). The Gibbsian factor of the potential U is the map e ^ : P([to]) —>• 

Definition 7 The connected function or Ursell coefficient (see e.g. US, I2f and references therein) 
is a map fip : P*([to]) -> I : I 4 cj)p(X), where cj>p(X), for all non-empty X C [to], is defined 
(uniquely and recursively) by the equations 

e -pU{X) = JJ 4>p(a) (3.2) 

Trgn(.Y) aGn 

Lemma 1 The following identities hold: 

M*) = e n ( e ~ mi -!) ( 3 - 3 ) 

g€G x {i,j}£E g 

where Gx denote the set of all connected graphs with vertex set X and for g e Gx, E g is the set 
of edges of g. 

Proof. Recalling (3.1), it is possible to invert the set of equations (3.2), i.e. to write down explicitly 
the Ursell coefficients. This is achieved by a so called Mayer expansion on the Gibbsian factor. In 
fact, for all finite X C N we can write 

e -m x ) = [(e-Wi _ 1 ) + 1 ] = Y, n M pVij - 1 ) 

{i,j}cx g&Gx {i,j}eE g 

where Qx denote the set of all graphs (connected or not connected) with vertex set X. Collecting 
together the connected components in the sum g x we S e t 

e -m x ) = y U K M) ( 3 - 4 ) 

7rgII(A') oGtt 

with 

Kp{a) = e n (f pVij - : ) ( 3 - 5 ) 

g&G a {i,j}£E g 

Comparing (3.2) with (3.4) we get Kp(X) = fffiX), for all X C [to]. □ 

Remark Note that fiXX) = 1 when |X| = 1 since the sum over Gx contains just one graph with 
no edges and the empty product is equal to one. 

Let Q n be the set of graphs (connected or not connected) with vertex set [to]. A graph g € Q n 
induces a partition n(g) € n([n]) whose blocks are the sets of vertices of the connected components 
of g. 
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We denote by T n the set of all trees with vertex set [n] (connected graphs in [n] with no loops, or 
also, connected graphs in [n] with exactly n — 1 edges). For a fixed tree r E T n with edge set E T , 
we denote by $ T the set of all bijections ip T : [n — 1] —>• E r (the set of all labelings of the edges 
of r). Given a tree r € T n and an edge labeling ip T E <F r we denote by the graph with vertex 
set [n] and edge set E Tk = {{i,j} E E r : ^ 1 ({i,j}) < k} (in other words t*, is the graph in Q n 
obtained from r by erasing all edges with labels greater than k). 

Theorem 3 (Basuev 1979) The following identity holds 


<MM) = 


e n 

g£G n {i,j}£E g 


-PVij 


- 1 


rP rPi r 

= (-ir“ 1 / ^ i/ dfh ■ ■ 

Jo Jo Jo 


dPn-l y II Vi 

r&T n {i,j}eE 


^ Yh e 

ip T e$ T 


E (Pk Ea67r(T fe ) ^( a ) 

k =1 


where /3 n = 0 by convention. 


(3.6) 


Remark . In the case of pair potentials with a hard-core (e.g. Vij = +oo for some {i,j}), one 
has just to introduce a cut-off H > 0 to give meaning to the expression in the r.h.s. of (3.6). 
One defines an auxiliary potential Vij(H) such that Vij(H) = V tJ if V. L j < +oo and Vij(H) = H if 
Vij = Too. Both sides of the identity (3.6) are well defined as long as the pair potential V(H) is 
considered. One then can just take the limit H —>• oo to recover the identity (3.6) for the hard-core 
potential. 

In the original paper [3] Basuev did not explicitly state the algebraic identity (3.6). For this reason, 
and also to make this paper as self-contained as possible, we present in Appendix A a detailed proof 
of the identity (3.6). Of course, all ideas contained in the proof are due to Basuev. 


4 Estimates for the Mayer coefficients 


In this section we show how to use the Basuev tree graph identity (3.6) introduced in the previous 
section to get direct and “combinatorially good” (in the sense of (2.9)) upper bounds on the Ursell 
coefficients (2.7). We will see in Section [5] that such bounds are better than all known previous 
bounds (from the old Penrose-Ruelle bound (2.18) to the recent Morais-Procacci-Scoppola bound 
(2.20)) as long as we restrict ourselves to the class of stable and tempered potentials introduced in 
Definition [2] which, we recall, embraces nearly all relevant examples of pair potentials in physics. 
As previously said, a Basuev pair potential according to Definition [2] has a property which is crucial 
in order to obtain these bounds. This property is expressed by the following theorem. 


Theorem 4 (Basuev 1978) Let V be a Basuev potential according to Definition [3J then V is 
stable. Moreover the representation 


with 


V a (r) 


V = V a + K a 

(V(a) 

if r < a 

1 V(r) 

if r > a 
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(4.1) 




and 


(V{r) — V(a) if r < a 

K a (r) = \ (4.2) 

[ 0 if r > a 

is such that the potential V a defined in (4-1) is also stable and it has the same stability constants 
B and B of the full potential V. 

Once again, to make the present paper as self-contained as possible, we give in Appendix B the 
original proof of Theorem 0] due to Basuev. 

This feature of a Basuev potential V (i.e. V = K a + V a with K a positive compactly supported 
and V a stable, tempered and absolutely integrable, with stability constants equal to the ones of V) 
matches very well with the simple and handy structure of the exponential factor exp{— Xsk=i (Pk ~ 
Pk+i)Ylaeir(T k )U((x)} appearing in the Basuev tree graph identity (3.6). It allows to use the 
condition of stability (2.10) to bound the part in the exponent coming from the term V a ■ At 
the same time, the positive term K a in the exponent controls the eventual divergence of the full 
potential at short distances. 

Indeed, from Theorems [3] and 0] the following estimate for the Ursell coefficients can be obtained. 
Once again, the following theorem is a development of a result due to Basuev [0]. 


Theorem 5 (Basuev 1979) Let V be a Basuev pair potential with stability constant B. Letn € N 
such that n> 2. Let C n (/3,A) be defined as in (2.7) and let B be defined as in (2.14). Then there 
exists a > 0 such that the following upper bounds hold 


1) 


C n (P, A) < 


n 


n—2 


n! 


2 n/3 B [C*(P)] n ~ l 


(4.3) 


with 


r (l — e -P[V(\x\)-V(a)\\ r 

C*(f3) = / dxP\V(\x\)\ \ > + / dx/3\V(\x\)\ (4.4) 

J\x\<a P[V(\x\) - V(a)\ ,/ N > a 


2 ) 


C n (P, A) < 


n 


n—2 


n\ 


e P B _ 1 

71—1 

r 



C(p,B) 

pB 


L J 


n—1 


(4.5) 


with 


_ r BB(1 — e —hW{\ x \)—V{a)—B]\ r 

C{(3, B) = / dx(i\V(\x\)\ J \ ' J + / dxfi\V{\x\)\ (4.6) 

J\x\<a j3(V{\x\) — V(a) — B)(ed B — 1) J\x\>a 


Therefore the series (2.6) converges absolutely, uniformly in A for any complex A inside the disk 


|A| < max 


1 

e P B+1 C*(/3) 


pB 


e(eP B - 1 )C(P,B) 


(4.7) 
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Remark. The first inequality (4.3) can be deduced also by using the early version of the Brydges- 
Federbush tree graph identity (e.g., the version presented in Theorem 3.1 in [6] or in Theorem 2 in 
m)- On the other hand, as said at the beginning of Section [3j the later versions of the Brydges- 
Federbush tree graph identity (such as the one given in Theorem VIII.3 in [8]) could have been 
used to obtain the second inequality (4.5). 


Proof. We first prove 1) i.e. inequality (4.3). Let us denote shortly 

^( X1 ,...,x n )= n 


g&Gn { i,j}£E g 

Starting from the Basuev tree graph identity (3.6) with V t j = V (|Xi — Xj |) 
the following upper bound 


we have immediately 


\(f)p(xi,...,X n )\ < 


r& rfi i r 

dp i / dp 2 --- 

Jo Jo Jo 


dPn-1^2 \V(\Xi-Xj\)\ 

T£T n {i,j}£E T 


X 


>V V] {Pk fik+l) 5Z Q! g 7 r('r fc ) U (a) 

Y e k=1 


VVG# T 


(4.8) 


By hypothesis, since V is Basuev, we have that V = V a + K a where V a is stable with the same 
stability constant of V and K a is positive and supported in [0,a]. Therefore, for any a C [n] we 
have that 

U(a)= V(\ Xi - Xj \) = U a (a) + U a (a) 

{i,j}Ca 

with 

U a (a) = V a(\ x i~ x j\) 

{i,j}Ca 

and 

U a (a) = Y K a (\xi-Xj \) 

{i,j} Ca 


Now, since V a is stable with stability constant B we have that, for all a C [n] 


Y V a (\xi-xj\) >-B\a 

{i,j}Ca 


Thus, for any partition ir of [n] 

U a (a ) > — B |or| > —Bn 

aG7r aG7r 

So we get (recall that, for all k € [n — 1] we have Pk ~ Pk+\ > 0 an d Y^k=\_(Pk — Pk+i) = Pi < P) 


n—1 


71—1 


exp i Y^ k - ^ fc+i ) Y r = exp {- Y^ k - P k + i ^ Y i u *( a ) + u a (a)\ > < 


k =1 


aG7r(r i .) 


k =1 


a£tr(T k ) 
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n— 1 

< exp { - J2(Pk ~ Pk+i )[~Bn + ^ U a (a)\ J> < e pBn e~ ^Zi(Pk-Pk + i) E Q6 „ (Tfc) u a {a) 

k=l ae7r ( Tfc ) 


We have now 


n— 1 


n— 1 


X>*-&+i) E Ua(a) =Ea[ E ^a(a)- X] ^a(“) 

aen(Tk) fc =l a€7r(rfc) ae7r(Tfc_i) 


fc=l 


The factor Eae^) ^(a) - Eoeih.O U a (a) is a sum of positive terms of the form if a (|x s — x;|) 
indexed by edges {s, Z} € E n including the edge {i, j\ of r which has the label k, that is, such that 
(f T {{hj}) = k. Therefore we get 


n— 1 


EC&-&+ 1 ) E ^( a ) - E ^diM’i)^ 

«e7r(r fc ) {2,j}e-Er 


a VI 


(4.9) 


k=\ 


Hence 


n —1 


ex P -E^-ft+i) U(a)\<eP Bn JJ e -^r({ij}) K a(\ x i- x j\) 


k =1 


o£Tt(r k ) 


Therefore we get 


r/3 r/h 


{i,i}G£; T 


r/3n-2 


< e nBP / dp 1 dp 2 ... dp n -i F((3i,...,(3 n -i) 

Jo Jo Jo 


with 


nsno = E E n i^(i £ 

reT n (/? T e^ T {2,j}e^r 


q A£ T ({£,j})-^-a(l a '£ ^jl) 


(4.10) 


(4.11) 


The sum over labels ip T makes the function F(/3 1 ,... ,P n — 1 ) in the r.h.s. of (4.11) symmetric in the 
variables Pi,... ,Pn— i- This means that in the r.h.s. of (4.10) every variable Pi can be integrated 
from 0 to P as long as one divides by (n — 1)!. Hence we get 


\<j>p{xi,...,x n )\ < 


0 nB/3 


rP 


,, ,,,E E n / 4 ^<{«)> 

1 B TGTnfre'Pr {i,j} & E T J0 


s nB/3 


(n — 1)! 


E n l F (l : 

rST„ {i,j}£E T 


-- p nBfi 


= e n 


Xi - Xj I 


tST„ {i,j}eE T 
i *t- 

jflA'adSi— 


(1 _ e -/3^o(|xi-xj|)\ 

jl)l ifa(|Xi-Sil) 

[1 _ e -P K a{\xi-Xj\)^ 

pK a (\xi - xj |) 


E i 




where of course 1 — gK f | ^ 1 = 1 when K a (\Xi — Xj\) = 0. Therefore 


1 1 


nB/3 


\C n (P, A) | < JTT 7 / dx 1 ... / dx n \(j>p(xi,...,X n )\< -j- ^ 

l A l n - J A J A n ■ — r-'-p 


W T 


r£T n 
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with 


w T = 


w\L dxi -I , dx \.n |v(| 


Xi - Xj I 


[1 _ e -PKa{\xi-Xj\M 

K a (\xi - Xj|) 


{*,j}e£V 

We can estimate w T , for any r E T n , as follows. Let us set, for e = {i,j} E E T , y e = — x j-, then 

I n —1 


UV < 


TTT / <^1 IT / ^ 

I a I7a 


fl _ e ~P K 4\ye\)] 
| [ *.(1*1) 1 


[ dx \V(\x\)\ 

JR d 


[1 _ e -/9#a(M)] 
A 'a(M) 


Therefore we get 


n —1 

£» - 

T£T n 


e nB/3 n n—2 

n ! 

where in the last line we have used the Cayley formula GT n 1 = n n 2 [9]. Finally, recalling 
definitions (4.2) and (4.4), the bound (4.3) follows. 

We conclude by proving 2), i.e. inequality (4.5). 

If we use the constant B we have that, for any a C [n] 

y, V a {\ x i - Xj\) > -B(\a\ - 1) 

{i,j}Ca 

Thus, for the partition ^ r E Ll(rfc) of [n], we have that 

E U a (a)>-B y (H-l) 

oG7r(r fc ) aGir(T k ) 


dx /3|W(|a:|)| 


[1 _ e -/3A-a(M)] 
/^a(M) 


|Cn(/3, A) | < 


0 nB(3 


n\ 


7 

J R c 


dx |C(|x|)| 


[1 _ e -/3ATa(|a;| 

^a(M) 


Since Tk has exactly n — k connected components, we have 

y (M - 1 ) = k 

CKG7r(-r fc ) 


So we get 

y U a (a) > —~Bk 


«G7r(rfc) 

Therefore, recalling that U{a) = U a (a) + U a (oi ) and that /3k — /3k+\ > 0 for all k E [n — 1], we get 
the following bound for the exponent in the integrand of the r. h. s. of inequality (4.8) 


exp < 

n—1 

-y(A-A + i) y u{a) 

> < exp < 


„ fc=1 ae-7T(r k ) 



n —1 


yc B k -p k + 1 )[-Bk+ y U a (a)} 

a£ir(r k ) 


k =1 
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Now observe that 


Pk+l)k — ^ ' Pk 


and, as before (see (4.9)) 


^(Pk-Pk+l) Y ^ Y Pvr{{i,j}) K ^\ x i- x i\) 

k= 1 o€7r(r fc ) {ij'}eS T 


so that 


71—1 71—1 _ _ 

- E (Pk-Pk+l) E U(a) E A®- E ^¥>t({*,j}) K “(I x *- x J'I) “ E /Vr({i,j})[ K “(l x «- x :*l)- s ] 


aeir(T k ) < gfc=l {ij}6£r 


— g 


Therefore plugging inequality above into (4.8) we get the estimate 


rh rPi rPn -2 ^ ^ 

| < /> /3 (.Ti,...,x n )| < / dpi dp 2 ... dp n -iY Y n 1^(1 

40 JO Jo ^ r, 


Xi-Xj\) I X 


rST„ (^T-G<I> r {*d}eE r 


E /^tCOo}) Ka(\%i Xj |) 7? 

X g {*o'}e^r L 

As before, we can symmetrize in the r.h.s. the integration domain of the variables Pi, ..., P n ~i 
integrating all them between 0 and P and dividing by (n — 1)!. We thus obtain 


\<t>p(xi ,..., x n )\ = Y II d ; 

rST„ {i,j}eE T 


[1 _ e -P(Ka(\Xi-Xj\)-B )j 

(Jt a (|xj - xj\) - B) 


eP B - 1 

PB 


E n 

rGT n {i,j}G^r 


/35[1 — 

'(^-^(ITadSi-SjD-B) 


where the factor 




^(ifadxj-XjD-S) 


—. 1 is set equal to 1 when K a (\xi — xj\) — B = 0. Therefore, 


r r~r 11 n ~i i 

N-] nn " 2 r f 1"“ 


!<?„(/?, A)| < 


[ dx P\V(\x\)\-^ 

jRd (e PB 


(e^ - l)P(K a (\x\) - B) 


Then, recalling definitions (4.2) and (4.6), the bound (4.5) follows. □ 

We stress that, in most of the cases, inequality (4.5) is stronger than (4.3). Let Pq is the solution 
of the equation (exp{/31?} — 1 )/PB = exp{/3L>}. Then the fact that (4.5) is stronger than (4.3) for 
all temperatures P < Pq follows from the following proposition. 

Proposition 1 For all \x\ < a and for all PB > 0 we have that 


pB(l — g-£[ y (M)-W“)]- B ) (1 — e -P\ v (\ x \)- v ( a )\) 

P(V(\x\)-V(a)-B)(eP B -1) ~ P[V(\x\) - V{a)\ 


(4.12) 
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Proof. Let us denote shortly 


(4.13) 


A = (3<y(\x\)-V(a)) 

By assumption we have that A > 0 for all |x| < a. Define the function / : 


with 


( y(ey A -l) if / o A 
(y-A){ e v- 1 ) 


f(y) = < 


1—e 


A 
A 

K e^l 


if y = o 
if y = A 


Let us check that f(y) is, for all A > 0, monotone decreasing in M. Since f(y) > 0 for all 
we have that 

dlnf{y) f'{y) 


dy f(y) 


< 0 => f'(y) < 0 


But 


dhi /(j/) _ d o y ^ _ i) _ _ j n ^ e v _ i)l = 

dy ay 

= y(y — A) — g(y) 


where 


g(y) = 


if d eM\{0} 


if 2/ = 0 

Thus the r.h.s. of (4.15) is negative if g(y) is increasing. We have 

1 e y 

g\y) = -o - 


y 2 (e y — l) 2 


so that 


g'{y) > o 


( e w _ 1)2 e 2 y _ 2e y + 1 

y z < - -=---= e y + e y - 2 


e y 

I.e. g(y) is increasing (and hence f(y) is decreasing) if 


y 


2 <e y + e~ y -2 = 2 


e y - e~ y 


- 1 


= 2(cosh(y) — 1) = 2 


y^_ / 

2 4! 


which is true for any y. 

Thus the function f(y) is decreasing and therefore, since /3B > 0, we have 

fWB) < /(0) 

which, recalling definitions (4.13) and (4.14), coincides with the inequality (4.12). □ 


(4.14) 

y e R, 


(4.15) 
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5 New bounds for Lennard-Jones and Lennard-Jones type poten¬ 
tials. Comparison with previous bounds 

In this final section we perform a comparison between bounds (4.3) and (4.5) deduced from the 
Basuev Tree-graph identity and the recent bounds for Lennard-Jones type potentials given in [24] 
and [T9] . 

Let us first introduce the standard definition of a Lennard-Jones type potential (a.k.a. inverse 
power low potential). See e.g. [13] . or definition 3.2.10 in [33] or formula (4.2.2) in [15] . 

Definition 8 A pair potential V on M. d is of Lennard-Jones type if there exist positive constants 
w, r\, r 2 , with r\ < r 2 , C, C', e such that 


V (r) > 


c 

y.d-\-e 

if r < r i 

—w 

ifr 1 <r<r 2 

C' 

if r >r 2 

y.d+e 


(5.1) 


A Lennard-Jones type potential according to the definition above is clearly tempered and it is long 
known that it is also stable (see e.g. [lanuimEaE])- Here below we prove that a Lennard-Jones 
type potential is also Basuev according to Definition (3[ 

Proposition 2 Let V be a pair potential on M. d such that there exist constants w, r\, r 2 , with 
ri < r 2 , and non-negative monotone decreasing functions f, g with domain in [0, oo) such that 


with 


and 


Then V is Basuev. 


' £(0 


V (r) > < 


—w 

—rj(r) 


if r < r i 
if r\ < r < r 2 
if r > r 2 


lim i(a)a d = +oo 

a —^0 


'\x\>X2 


dx r)(\x\) < Too 


(5.2) 


(5.3) 


(5.4) 


Proof. By Theorem 0] we just need to show that there exists a such that (2.15) and (2.16) are 
satisfied. Fix a € (0, ri), let w = max{u>, r]{r 2 )} 


r/(r) 


r/(r) 

if r > r 2 

id 

if r < r 2 
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Then, by construction fj is monotone decreasing and such that f Rd dx fj( \x\) < oo. Moreover by 
conditions (5.2) we have that 

V~ < fj 

Hence, recalling (2.17) and considering also that, since we took a € (0, n), by hypothesis H“(|x|) = 
0 for all \x\ < a, we have 

n 

fi(a) < sup „ E^i) 

nGN, (xi,...,xn)eR dn ^^ 

\x{-Xj\>a, | Xi\>a 

To bound from above X^=i '/(I x i\)i having in mind that all particles are at mutual distances greater 
than a and at distance greater than a from the origin, we proceed as follows. We draw for each 
Xj a hypercube Qj with side-length a/2\fd (such that the maximal diagonal of Qj is a/2) in such 
a way that Xj is a vertex of the cube Qj and at the same time is the point of Qj farthest away 
from the origin 0. Since any two points among x \,..., x n are at mutual distances > a the cubes so 
constructed do not overlap. Furthermore, using the fact that fj is monotone decreasing, we have 

v(\ x j\) < f dx vi\x\) 

a JQj 


Recall in fact that the cube Qj is chosen in such way that \x\ < \xj\ for all x € Qj. Therefore 


i= 1 



dx r/(|x|) 




dx r/(\x \) < 



C d 


a 


d 


where C d = ^r~ f K d dx r](\x\). Hence we get 

, w c d 

^(°) ^ 

In view of condition (5.3), we can always choose a such that £( a)a d > 2 C d . Thus we get 

H(a)a d > 2 C d e(a) > 2^ =► V{a) >2^ => V(a) > 2 fi(a) 

a a a a 

□ 


Corollary 1 A Lennard- Jones type potential according to Definition^ is Basuev 

Proof. Just observe that setting £(|x|) = C/\x\ d+e and rj( |x|) = C'/\x\ d+e a Lennard-Jones type 
potential satisfies the conditions (5.2)-(5.4). □ 

In what follows we first compare the bound given by (4.3) with the Morais-Procacci-Scoppola 
bound (2.22) and show that bound (4.3) always beats bound (2.22). As we have seen, if B = B 
or the temperature is not too low, bound (4.5) beats (4.3). We then consider the specific case of 
the classical rescaled Lennard-Jones type potential V(r) = ~ri ~ /f- In this specific case, due to 
known computational results it is known that B and B are, if not equal, then very close to each 
other (i.e. B < (1.001)H). So, in order to compare our results in this case with those obtained by 
de Lima and Procacci in m we set, accordingly to them, ft = 1 and we show that (4.5) strongly 
beats the de Lima and Procacci bound. 
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5.1 Lennard-Jones type potentials 

As we have shown in Corollary |U a Lennard-Jones type potential V according to definition [8] is 
Basuev. I.e., if B is the stability constant of V, then there exists a > 0 such that V = V a + I\ a 
with V a stable and absolutely summable with the same stability constant B of the full V and I\ a 
positive supported in [0, a] (V a and K a defined as in (4.1) and (4.2)). 

Let us thus use this information to compute the (lower) bound given by Morais et al. for the 
convergence radius, say R, of the Mayer series of a gas whose particles interact through a Lennard- 
Jones type potential V. According to Theorem [2] this bound is as follows. 


eP B+1 C(/3) 


(5.5) 


where 


cm 



1 _ e -fi[V(\x\)-V(a)] + £y( 0 ) 



dx /3|V _ (|x|)| 


(5.6) 


Note that, in force of Corollary [1} we have replaced in the r.h.s. of (5.5) the (generally larger) 
stability constant B of the absolutely summable potential V a (appearing in the original Morais et 
al. bound (2.22)) with the (generally smaller) stability constant B of the full potential V. 

Let us compare this bound with the bound, say R*, on the same radius derived from (4.3) obtained 
in Theorem 0 which is 


e P B+l C*([3) 


(5.7) 


with 


c*m = 


f (\ — e -P[V{\ *|)-V(a)h r 

/ dx/3V{\x\) { — , ,, + 

l\x\<a P[V(\x |)-V(a)J J\ x \> a 


dx 




(1 - e -W(\x\)~ v (a)}) + p V (a) 


So we get 


Cm-C*iP)=PV(a) / dx 

J |a;|<a 


dx f5 |V'(|x|)| = 

(l - e -PlV{\x \)-%)])' 
(5[V(\x\)-V(a)] 

1 _ e -p[V(\x\)-V{a)\ 

P[V(\x\)-V(a)] 


+ 


1 - 


|^|^ a 


> 0 


dx /3|C(|x|)| 


Indeed, recalling that V(r) > V(a) for all r < a and setting z = /3[V(|a:|) — 14(a)], the function 
1 — — z — is greater than zero for all z > 0. Hence we always have that 

R* > R 


I.e. the bound (4.3) is always better than the bound (2.22). 

Let us briefly consider also the case a pair potential with an hard-core, i.e. a potentials V/,. c (r) such 
that Vhc(r) = Too for r < a (with a > 0). As mentioned in Section [2l a tempered pair potential 
with an hard-core is Basuev, and, as mentioned in Section [3l potentials with a hard-core can be 
treated with the Basuev tree graph identity (3.6) by introducing an auxiliary potential V(H). For 
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a potential I4 C with hard-core as above the factors C*(/3 ) and C(f3,B) appearing in bounds (4.3) 
and (4.5) become 

C*(J3) = W a (d ) + 0 [ dx \V hc (\x\)\ (5.8) 

J |rr|>a 

C{P,B) = -^—W a {d) + P [ dx\V hc (\x\)\ (5.9) 

eP B - 1 J\x\>a 

where W a (d) the volume of the sphere of radius a in d dimensions. Inserting (5.8) into (4.3), one 
gets the same lower bound for convergence radius obtained in [29, [21] for potentials with a hard¬ 
core (see e.g. Theorem 3 in [21]). It is also clear that (5.9) inserted into (4.5) yields an improved 
lower bound for the convergence radius respect to that of Theorem 3 in [21]. 


5.2 Lennard-Jones potential 

In this and the next subsection we consider the specific case of the three-dimensional classical 
(rescaled) Lennard-Jones potential: 

y ( r ) = ^l2-^ ( 5 ' 10 ) 

We denote by B hJ its stability constant. 

We first compare the bound stemming from Theorem [5] with the bound given in [19] (see Theorem 
3). To be coherent with [19] we set the inverse temperature at the fixed value /3 = 1. 

In [T9] the authors show that, by choosing a = 0.3637, the potential V defined in (5.10) can be 
written as V = V a + K a with V a (r), defined as in (4.1), stable with the same stability constant 
B lj of the full Lennard-Jones potential V and K a (r ) = V{r) — V(a ) for r € (0, o] and K a (r) = 0 
elsewhere. The lower bound obtained in m, for the convergence radius, say -RlPj was deduced 
from Theorem [2] of Section [2] Namely, from inequality (2.22) with /3 = 1 one gets 


-Rlp 


1 

C'(l)e B Lj+ 1 


For the purpose of comparison with Basuev bound, observe that 


(7(1)= [ dx \\e~ KM — 1| + |K(|5 

J R 3 L 


> 


= 4-7T 


L 


0.3637 


dr 


1 


|0.3637| 12 |0.3637| 6 

> 37444 + 12381 > 49825 


dx |14(M)| = 

/ 


r 2 + 47t 


dr 


1.3637 


> 


_ -(b lj +i) 

Therefore the De Lima-Procacci lower bound 7 ?lp obtained in |19] can be at best e 


49825 > 


i.e. 


Rlp < 


e ( b lj +1 ) 


49825 


(5.11) 


Now we use the estimate (4.5) (which in this case is tighter than (4.3)) to obtain an alternative 
lower bound of the same convergence radius. Since we want to use (4.5), to compare efficiently 
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this bound with the previous one (5.11) we need to estimate £> LJ , which appears in (4.5), in terms 
of L> lj . To do this we take advantage here of the considerable amount of computational rigorous 
results on optimal Lennard-Jones clusters in the literature. In particular the values of the global 
minima for configuration with n particles has been calculated up to n = 1610 (see e.g. [39] and 
reference therein). The tables in [39] show that B n = -^jB n (we recall that B n , according to 
(2.12), is the absolute value of the minimal energy of a configuration with n particles of a Lennard- 
Jones gas divided by n) is less than i?iooi for all n < 1000. So using these data one can conclude 
that B hJ = sup n B n < (1.001).B LI where recall that B hJ denotes the stability constant of the 
Lennard-Jones potential (5.10). 

Therefore, setting /3 = 1 and a = 0.3637 coherently with above, the lower bound, say R, of the 
same convergence radius is, according to estimate (4.5), given by 


R = 


(1-001 )B h 


eieiroKj-li C(1,B] 


= h(B L 


e 


(5.12) 


with 


and 


h(u) = 


(1.001)u 


e iooo — e 


for u > 0 


(5.13) 


f B (1 - e -l v (.\ x \)-V{a)-B h:j }} 

/ dx nw) j , + 


* |ai|<a 


{V(\x\)-V{a)-B w )(e B m - 1) J\ x \> a 


dx 114(1x1)1 


First observe that the function h(u) dehned in (5.13) is increasing up to a value around u = 1000 
and it is known that 8.61 < B hJ < 14.316. The lower bound was first obtained in m (see also 
Lemma 3 in [36]) while the upper bound has been recently tightened by Yuhjtman [38] (the previous 
best upper bound [36] was B hJ < 41.66). Hence h(u) is increasing in the interval [8.61,14.316] and 


so h(B hJ ) > h( 8.61) > 8.69. Therefore we get from (4.5) that R is at worst -g — 


8.69 


e S LJ + 1 C(l,R 


i.e. 


R > 


8.69 


e S Lj+ 1 C'(l, H lj ) 


By bound (4.12) in Proposition |T] we have 


C(1,H lj )<C(1, j B lj =0)= / dx V(\x\) 
Computing explicitly C(1,H LJ = 0), we get 


(l — e -[V(.\ x \)~V(a)]\ r 

1 v( , h v(] ’ + dx\V(\x\)\ 

H(|x|) - H(a) J | X |> 0 


C(1,H LJ =0) = 4vr 


l 


0.3637 


dx 


Ixl 10 |x| 4 


+ 47T 


l 


dx 


0.3637 


Ixl 10 |x| 4 


l—g L|a:| iJ |a;| b 10.36371^^ lO-SeSTI 1 

I 2 I , 2 

frp I7F |0.3637| 12 “ l " |0.3637| B 

< 0.823 + 12381.1 < 12382 


+ 


(5.14) 
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Hence we have that 


R > 


8.69 


12382 • 


Bt.t + 1 - 


O (RlJ+ 4 ) 

1425 


so that for the ratio Rlp /R we have 


-RlP 

R 


l 

49825-e a Lj +1 

1 

885-e s Lj +1 


1425 ^ 1 

49825 ^ 34,96 


(5.15) 


In conclusion we obtain that the bound (4.5) produces a lower bound for the convergence radius 
which is, for the classical Lennard-Jones potential, nearly 35 times better than the de Lima-Procacci 
lower bound for the same convergence radius. 


5.3 Optimal bound for Lennard-Jones potential 

The bound (5.15) can be strongly improved by trying to find the optimal a which minimizes the 
factor C (which depends on a). This optimal a coincides with the a as large as possible. This is 
because the first of the two terms in the r.h.s. of (5.14) increases very slightly as a increases, while 
the second term in the r.h.s. of (5.14) decreases very fast as a increases. The total effect being that 
C(f3, B hJ ) decreases rapidly as a increases. Of course we cannot increase a as we like. According to 
Theorem HI a is constrained to satisfy the inequality 2fi{a) < V(a). Any a satisfying this inequality 
will do and, due to the previous discussion, the larger is a, the better. To find this optimal a, we 
need an upper bound as tight as possible for /u(a). For this purpose we can use the bound recently 
obtained by Yuhjtman (see [[38], Proposition 3.1 II). 


Proposition 3 (Yuhjtman 2015) Let V be the Lennard-Jones potential defined in (5.10), let a 
be such that 0.6 < a < 0.7 and let p,(a) be defined as in (2.17). Then the following bound holds. 


p{a) < 


24.05 


Thus the inequality 2 fi(a) <V(a) = ~ri — Jr is satisfied if there exists a E [0.6,0.7] such that 


24.05 1 

2 -^r < — 

aP a 1 - 


which is true as soon as a < 0.6397. Redoing the calculations in (5.14) replacing a = 0.3637 with 
this new optimized value a = 0.6397 we get 


C(1,R lj = 0) = 4vr 


L 


0.6397 


dx 


1 


M 10 |x| 4 


+ 47T 


roc 

/ dx 
J 0.6397 


I rj. I 10 I X I ^ 


_ g L | n :| 12 | x | t > | 0 . 6397| 12 | 0 . 6397| 6 J ) 

~n 2 i i 2 n 

ijxp \x\f I0.6397P T 10.63971° > 

< 2.5 + 61.63 < 64.13 


+ 


With this bound for C we get our optimal lower bound for the convergence radius of the Mayer 
series of the Lennard-Jones gas at = 1 which is 


- 8.69 e p_ ( b lj +1 ) 

R > - > 


64.13 


7.4 


(5.16) 
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which is greater than the previous best lower bound (5.11) given in [T9j by factor 6, 7 x 10 4 . 

Finally, it is worth to remark that, using the recent upper bound for the stability constant of the 
Lennard-Jones potential B hJ < 14.316 obtained by one of us [38], the lower bound of convergence 
radius of the Mayer series of the Lennard-Jones gas given in (5.16) is, in absolute terms, 

6, 7 x 10 4 x e 41 ' 66 " 14 ' 316 > 5 x 10 16 

times better than the previous known best bound available in the literature, i.e. the bound (5.11) 
given in [R] where the upper bound B hJ < 41.66 was used. 

Appendix A. Proof of Theorem \3\ 

We begin by obtaining an expression of the Ursell coefficients (fp{X) (see Section [3j Definition [7] 
formula (3.2)) which is alternative to the standard one given in LemmalU equation (3.3). 

Lemma 2 The following identities hold for all non-empty X C [n]. 

4>p{x)= (-i) |7r|_1 (M - l)! e -^^R«) (A.l) 

Tre npo 

Proof. The identity (I A. ID can be deduced directly via the Mobius inversion formula. We first 
recall that the set of partitions IJ(Jf) of a finite set X is a lattice, where the partial order < is the 
usual partition refinement. Namely, if tt, n' € n(X) then we say that n is a refinement of n' and we 
write n < tt' if any block a of vr is such that a C a! for some block a! of tt' . Hence the maximum in 
the lattice n([n]) is the one-block partition id = {[n]} while the minimum is the n-block partition 
@ = {{l},{2},...{n}}. 

Let / : n(X) -> K be defined such that 


/(tt) = II ^( a ) 

aE7r 


Let us further define g : n(A") — > R such that, 

g(n) = [] 

QE7T 


Then, recalling that id G n(X) is the one-block (maximal) partition of X, g(id) = e hU{x) an( ^ 
f(id) = <fp(X). So with these notations (3.2) rewrites to 

e-mx = m = ^ /(7r) 

7r< , $ 


Therefore, by Mobius inversion formula, (see e.g. m) 

(/>p{x) = f(id) = hix-, v)g{K) 

7r<i9 
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where n(n, i?) is the Mobius function of the lattice of set partitions. It is now well known that 
/j{i r, i?) = (—1)I 71 "! - 1 (| 7 r| — 1)! (see again [32] or also [3D], formula 5C.20), whence (IA.1D follows. □ 

Remark. Lemma [2] implies the non trivial identity 


mm) = y (- i ) ki i (i yr i—!) !e /3 e “ 6- c/(q ) = y n 

7 reII([n]) g&G-n {i,j}eE g 

In particular, if K is the number of partitions of the set [n] with exactly k blocks (i.e. = 

|n*([n])|), then by setting /3 = 0 in (1A.2I) . it follows that 

n 

£(- 1)‘ _1 (*-W = 0 (A.3) 

k=l 

The representation (IA.2D of the Ursell coefficients is the starting point in order to prove the tree 
graph identity given in 1979 by Basuev 0], Moving on towards this goal, let us introduce some 
further notation. We set shortly 

rP rPi rfin -2 r 

/ d(5\ I d/3 2 ■■■ dj3 n -1 = d(3 

Jo Jo Jo Js 

Let A = {ai,..., a*,} be a family of disjoint non-empty subsets of [?z]. An element a of the family 
A is called a block and k = |A| is the cardinality of the family. We let q(A) = {a C A : |oj = 2}. 
Namely q(A) is the set of all unordered pairs of blocks of A. Note that a E q(A) implies that 
<7 = {a, a '} with a C M, a> c M both nonempty and such that aC\a' = 0. Moreover, for a E q(A) 
and A family of disjoint of non-empty subsets of [n], we denote 

W a = Y Vi i ( A - 4 ) 

{ i,j}C[n] 
i€o:, j€ot f 

U(A)= Y W ° ( A -b) 

treq(A) 

Note that (IA.4I) and (IA.5I) jointly with Definition [5] imply that U(A) = D(U^ = 1 a i) — Y2i =l U(oti) 
and also U(a) + U(a!) + W a = U(a U a'). Note also that if |A| = 1 (i.e. if A contains just one 
block) then U(A) = 0, since q{A) is empty. 

Let 7 = {ai,..., afc} be a partition of [n] and let a = {a,, aj} E q( 7 ), then we denote by P a 7 the 
new partition of [n] given by P a 7 = (7 \ { oti,otj }) U {ai U aj}. Namely P a 7 E LI([n]) is the new 
partition of [n] given by 

Pa 7 — { 07 ,..., aj— 1 , a ^|_i, • • •, aj— 1 , aj+i ,..., aai U aj } 


o-PVi 


ij — \ 


(A.2) 


Note that \P a l\ = M — 1- 

Finally, if 7 E LI([n]) then we denote by 1 ^( 7 ) the set of all partitions of 7 with fixed cardinality l. 
We stress that 7 E LI([n]) is a partition of [n]. So an element n E 1 ^( 7 ), is in general not a partition 
of M but rather a “partition of a partition”. Namely, a block A E ir is in general a disjoint family 
of sets of M- 
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With the notations above, for any partition 7 E II([n]), we define 


M 

^> 3 ( 7 ) = ^^(— l) fc-1 (fc — 1)! e -PI2AeirU(A) (A. 6 ) 

fc=i 7ren fc (7) 


Note that 

#( 7 ) = 1 if l7l = 1 (A.7) 

Observe also that, by (1A.3D . we have, for any 7 E II([n]) 

0/3 ( 7 ) = 0, for (3 = 0 (A.8) 

Observe finally that, with these notations, if we denote by 70 the partition of [n] given by 
70 = {{!}, {2},... , {n}}, then by equation (IA.2D we have that 


^(to )== y n 

g£G n {i,j}&E g 



Lemma 3 For every 7 E IT([n]) with I 7 I >2 we have 

0 / 3 ( 7 ) = - f d$ Y W a e-^' w ^p(P al ) 

0 0 - 69 ( 7 ) 


(A.9) 


Proof. By definition (IA.6D we have that 0 ^( 7 ) is differentiable as a function of /3, hence, using 
also (1A.8I) we can write 



M 

Ef-iN 1 /*-!) 1 E 


Y U{A)e~P'^^ 


k=l 


TTGn fc (7) AG-it 


Hence we just need to show that 


ItI 

Y W c e- pw ^ p {P^) =^(-l) fe - 1 (jfc- 1)! Y Yu(A)e-^ A e^(A) (A.10) 

crEg(7) k=1 'irEllk('y) A(El7t 


We have 


H-! 

Y Wae-P^ippiPrf) = Y(- l) fc_ 1 (fc-l)! Y W ° e-^A^U(A)+ w ^} 

oeg( 7 ) k= 1 iren fc (P CT 7) 

Let 7 T be a partition of 7. We write a 1 Z 7 r if there is A E 7 r such that a C A. With this notation 
we can write 

Y e -h[T.A^U{A)+ w a ] = J2 e -^AenU(A) 

3 ren fc (P CT 7 ) ’ r en fc ( 7 ) 
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Hence 


Y W a Y c~Ma^U{A)+wA [ = 22 W ” 'll e ~^A^U{A) = 


0-69(7) ifGn fe (p cr 7) 


7 ren fe( 7 ) 

cr[l 7 r 


= ^ e -PTAMt A ) 

E ^ 


E^) 

^en fc ( 7 ) 

0-69(7) 

cr|Z 7 r 

7 T-en fe ( 7 ) 

AG 7 T 


Note that Yla & q { 7 )- aon W a = Xm6tt ^'(-' 4 ) because both sides are the sum of W a over all pairs of 
blocks in 7 such that both blocks are in the same block A in 7r. Hence we get 

ItI— 1 _ _ 

Y W' T e-P Wa xl) fi {P <r 7)= j](- 1 ) fc “ 1 (fe- 1 )! 22 e-PZA&rUW ^ U{A) 

a £q( 7 ) k=l TT^H] t ( / y) AGtt 

When k = I7I there is only one partition 7r in I7I blocks and 7r is such that each block A E 7r 
contains just one element and thus U(A) = 0 . Hence we can change the upper limit I7I — 1 in the 
sum over k to I7I in the r.h.s. of equation above so that it becomes ( 1 A. 10 D . □ 

Definition 9 Let k € [n], we define a class &k of sequences (07,..., Ok) as follows: (07,... , 07) € 
6 if <J\ E 5(70) and, for i > 2 , ai E q{P ai ^i ■ ■ ■ P ai 70) • Given (07 ,... ,07) E 6 we also denote, for 
all i € [k], 7 i = P r 7 i P(T i _ 1 ■ ■ ■ Pa i7o- We finally denote shortly & n -i = 6. 

Corollary 2 Let V be a pair interaction in [n]. Then the Ursell coefficient <f>p([n]) satisfies the 
following identity. 


MW) = MM = MM 1 [dp 22 w *i ■ ■ • 

S (<7i,...,<7 u _i)e6 


Proof. Just iterate n — 1 times the equation ( 1 A. 9 I) to fig(jo) and observe that by construction 
I PfT n -i ■ ■ • Pa i 7 o| = 1 , so by (lAJj), 


P/ 3 (Pan-1 ■'■Pax 70 ) = 1 

Lemma 4 Let k € [n] and let (07,..., 07) E 6^, then 

k 

22 m = 22 u m (A.11) 

i= 1 o€ 7 fc 

Proof. By induction on &. For k = 1 the identity ( 1 A. 11 D is true. Indeed, supposing 07 = {i, j} 
and hence 

7 i = {{!}, ••• ,{*-!},{* + !}:••• ,{J - !}>{ 3 + !}>•■ ■, W 

we have 

w ai = v ij = u({i,j}) = u({i,j})+ 22 um = 22 u M 

l€.[n]: CtG 71 
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Assuming thus ( 1 A. 11 D true for k, let us prove that (1 A. 1 If) also holds for k + 1 . Let <r fc+1 = { aa"} 
with a' E and a" E 77. Then we have 

k-\-1 k 

E w °i = E ^+ w «> + i = E +"Vi = 

«=1 i=l 


E 17(a) + C/(a') +^(a") + ^ fc+1 


E 17(a) + *7(a' U a") = E ^(«) 

Q'7^Q! / , OL^hOL 11 


□ 


Corollary 3 Let 07,... , <r n _i 6e such that a % E g(7j_i) and 7* = P CTi 7j_i. TTie following identity 
holds: 

n—1 n—1 

Ew. = E(a-a+i) E tf(“) ( A - 12 ) 

i=l fc=l «G7i! 

where we have put ( 3 n = 0 . 

Proof. We have that Pi = YlkZiiPk — Pk+ 1) and thus, using flA. 111 ) 

n—1 n—1 n—1 n—1 k n—1 

E = E E(& - a+i)^ = E(& - &+i) E ^ = E0 3 * - a+i) E ^ a ) 

i=l i=l fc=i fc=l i=l fc=l 067 k 


□ 

To state the following lemma we need to introduce some further notation. Recall that Q n is the set 
of all graphs with vertex set [n] and, given g E G n , ir(g) denotes the partition of [n] induced by g 
such that a block a € ir{g) is formed by the vertices of a connected component of g. 

Given 7 E II([n]) and a = {a, a'} E 5(7), and given A = {i,j} E E n , we write A < <r if i E a and 
j E a'. 

Lemma 5 Let T n be the set of all pairs ((07,... ,cr ra _i), (Ai,..., A n _i)) with (<ti, ..., <r n _i) E © 
and, for i = 1 ,..., n — 1 , A; E E n and A* <1 cr*. Let T n be the set of all pairs r = (r, cp T ) where r is 
a tree with vertex set [n] and (p T is a labeling of the edges of r (i.e. a bijection (p T : E T [n — 1 ] : 
{i, j} i->- <p T ({hj}) with E t denoting the edge set of t). Then there is a bijection between S n and 
T n . 


Proof. Given the pair ((07,... , <r n _i), (Ai,..., A n _i)) E T n , r = (Ai,..., A n _i) is a set of labeled 
edges in [n]. Let us prove that r is a tree with vertex set [n]. Consider the graph 17 = (Ai,..., \k) 
(with k < n — 1 ). By construction the connected components of 77 are the blocks forming 7^ = 
Pa k •••P o - 1 70 ) i- e - vr(rfe) = 77. Hence r = r n _ 1 has only one connected component, i. e. r is a 
connected graph with n — 1 edges, thus a tree. The labels of the edges are the subindexes, so 
T — (Ai, . . . , A n _i) E Tn- 

Conversely let r = (r, <^ r ) E T n with labeled edges (Ai,..., A n _i) and let 77 be the graph in Q n 
obtained from t by erasing all edges with labels greater than k. Then a unique sequence 71,..., 7 n _i 
of partitions of [n] is determined by posing 7*. = it (77). This sequence 7i,...,7n-i uniquely 
determines a sequence 07, ..., <r n _i such that 07 E q{p/i) and 7j = P ai Ti -1 and by construction 
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A i <\ (Tj . Hence we have constructed a function that associates to an edge-labeled tree r € T n 
with labeled edges (Ai,..., A n _i) a unique pair ((<ti,..., (Ai,..., A n _i)) € T n . We call 

(cfj",..., &n_i ) this unique sequence determined by the edge-labeled tree r = (Ai,..., A n _i). □ 

Proof of (3.6). By Corollary [ 2 ] we have 


<?V(M) = ^( 70 ) 



then ( 3 . 6 ) follows by Corollary 0 formula ( 1 A. 12 IK 


Appendix B. Proof of Theorem [4] 


We start recalling that 


V~(r) 


\V(r)\ - V{r) 


The thesis is trivial if V =0 (i.e. if V is purely repulsive). So we may assume that V 7^ 0 . For 
any (aq,..., x n ) € M nci and any i € [n], let 


Ei(x i,...,x n )= ^2 V(\xi-Xj |) 

je[n]:j^i 


so that 

U(x i, ...,x n ) = Ei(x i, ...,x n ) + U(x 2 , ...,x n ) 

Let now (x\,... ,x n ) G M nci be a configuration in which there is a particle, say in position x\ 
(without loss of generality), such that E\{x \,..., x n ) > 0 . Then we have 

U(x i, ...,x n )> U(x 2 , ...,x n ) 

therefore, since - < — 

’ n n —1 ’ 

-~U(xi, ...,x n )< - U(x 2 , ■ ■ ■ ,x n ) 

n n — 1 


Thus we have that the configuration (x2, ■ ■ ■ ,x n ) produce a value —U(x 2,... ,x n )/(n — 1 ) which 
is nearer to B than —U(xi,... ,x n )/n. Whence we can look for minimal energy configurations 
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(xi,..., x n ) limiting ourselves to those configurations in which the energy per particle Ei(x i,..., x n ) 
is negative for all i € [n]. 

Now let us consider the system of particles interacting via the pair potential V a defined in (4.1) 
and let us assume that conditions (2.15) and (2.16) holds. Note first that, due to condition (2.15), 
V~(\x\) = max{0, — I4(|x|)} = U _ (|x|). Consider then a configuration (xi,... ,x n ) such that there 
exists {i,j} C [n] such that |x, — Xj\ < a, thus there is at least a particle, (which, without loss 
of generality, we can assume to be the particle indexed by 1 at position xi = 0 ), which has the 
maximum number of particles among xz,...,x n at distance less than or equal to a. Say that the 
number of these particles close to xi less or equal to a is l (l > 1 by assumption). The energy E\ 
of the particle at position xi can thus be estimated as follows. 

Ei{x i,... ,x n ) > IV(a) - £ w(M) 

fc£[n] 

\ x k I > a 


To control the sum U~(|x|) observe that we are supposing that each particle has at most l other 
particles at distance less or equal than a. Thus take the k € [n] such that V~(\xh\) is maximum. 
Again, without loss of generality we can suppose k = 2. In the sphere with center X 2 and radius a 
there are at most l + 1 particles (the particle at position X 2 plus at most l other particles). Hence 

E V-(\x k \)<(l + l)V-(\x 2 \)+ E y '(NI) 

k£[n] k£[n] 

I > a \xk\>a,\x k -x 2 \>a 


Iterating we get 

E H-(|x fe |) <(/ + !) E ^“(W) (B- 1 ) 

\ x k\>a \ x i~ Xj\>a 

where in the sum in the r.h.s of (IB. ID all pairs of particles are at distance greater than a to each 
other. Therefore, recalling definition (2.17) we have that 

E B“(|x fc | )<//(o) 

ke[n] 

\ x i~ x j\> a 


and hence 

Ei(x i,..., x n ) > IV(a) - (l + l)fi(a) 

so we have E\ > 0 whenever 

v ( a ) > 

Using assumption (2.16) and since < 2 we get 

Ei(xi, ... ,x n ) > 0 

In conclusion, if a configuration (xi, ..., x n ) is such that some particles are at distance less or equal 
than a, then there is at least one particle whose energy is positive. Hence the minimal energy 
configurations for V a must be searched among those configurations in which all particles are at 
distance greater than a from each other. But for these configurations V a = V which implies that 
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V a and V, if stable, have the same stability constant B (and similarly also the same B ). Finally, 
it is easy to see that V is stable. Just observe that for any configuration (aq ,... ,x n ) for which 
particles are at distance greater than a from each other we have 

1 n 1 

U a (x 1 ,..., x n ) = 2^2 V( y\ x i~ x i\) > -T^Ao) 

i= 1 je[n] 


which implies that V a (and a fortiori V) are stable with stability constant B < Ay-. 
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